SU-4252-816,TKYNT-05-19, IMSc/2005/09/21 



Semi-superfluid strings in High Density QCD 

A. P. Balachandran 1 , S. Digal 2,3 , and T. Matsuura 2 
1 Physics Department, Syracuse University, Syracuse, New York 132^-1130, USA 
2 Department of Physics, University of Tokyo, Tokyo 113-0033, Japan 
3 Institute of Mathematical Sciences, C.I.T. Campus, Taramani, Chennai 600 113, India 

We show that topological superfluid strings/vortices with flux tubes exist in the color-flavor 
locked (CFL) phase of color superconductors. Using a Ginzburg-Landau free energy we find the 
configurations of these strings. These strings can form during the transition from the normal phase 
to the CFL phase at the core of very dense stars. We discuss an interesting scenario for a network 
of strings and its evolution at the core of dense stars. 
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I. INTRODUCTION 



Color superconductivity is expected to be the ground state for high baryon chemical potential 
/i and low temperatures 0. It is believed that such a state of matter may exist at the core of 
very dense stars. To find any signature of its existence, it is important to study various properties 
of the color superconducting phase. When all the quarks are masslcss, there are superconducting 
phases, namely the color-flavor locked (CFL) and the 2SC phasespj. In the CFL phase, all quark 
flavors take part in the condensation, but in the 2SC phase, one quark does not participate in the 
condensation. With realistic values of quark masses and charge neutrality conditions, the CFL phase 
gets modified to phases known as mCFL, dSC 0, uSC 0, gCFL, g2SC HQ, and FFLO phases 
[a, Li e tc- One of the interesting properties of a color superconductor is that it is also a superfluid. 
This is because apart from local symmetries, such as SU(3)c of color symmetry and U(\)em of 
electromagnetism, certain global symmetries are also spontaneously broken in the superconducting 
phase. The global symmetries spontaneously broken in the superconducting phase are SU(3)f (the 
flavor symmetry) and U(1)b (the baryon number symmetry). One expects defect solutions such as 
flux tubes or vortices to be present in the superconductor since it is both a superconductor and a 
superfluid. Recently there have been lots of studies of superfluid vortices 00 and flux tubes in 
the color flavor locked (CFL) phase and in the 2SC phase 11]. In these studies, superfluid vortices 
are topologically stable. The flux tubes studied so far are not topological and it is not clear if they 
are stable at all. 

All the topological string solutions like vortices or flux tubes are due to and are related to the 
existence of nontrivial loops (NNL's) in the order parameter space (OPS). An OPS is the set of 
all possible values of the order parameter (OP) which is the diquark condensate in this case. Each 
NNL in the OPS corresponds to nontrivial defect solutions in the superconducting phase. Previous 
studies have considered NNL's in the OPS which are generated only by the baryon number charge 
00. However when symmetry groups such as color, flavor and baryon number are spontaneously 
broken, one can have other NNL's. There are NNL's in the OPS which are generated partly by the 
baryon number and partly by other non-abelian color or flavor generators |l2j . These loops give rise 
to non-abelian strings. These string defects have been studied previously in the context of broken 
chiral symmetry in the framework of the linear sigma model 13]. In this work we consider the 
non-abelian defects in the color superconducting phase. We include the effects of gauge fields which 
result in flux tubes. These flux tubes are unlike those in ordinary superconductors. The energy per 
unit length of the string behaves like that of a superfluid vortex. This is why we call these non- 
abelian flux tubes as semi-superfluid strings. In this work we consider the semi-superfluid strings in 
the CFL phase. We will argue that these defects and superfluid vortices are possible in the dSC, uSC 
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and mCFL phases. In order that the flux tubes are dynamically stable, the color superconductor 
must be type II. For asymptotic values of the chemical potential /i, the color superconductor is type 
I, but for intermediate values of fi it can be type II (3, L1J] • 

This paper is organized as follows. In the next section, we discuss the Ginzburg-Landau free 
energy and discuss the nontrivial loops in the OPS which give rise to the non-abelian string defects. 
Section III will contain numerical results for the configuration of string defects in the CFL phase. 
In section IV we discuss a possible scenario for the string network and its evolution at the core of a 
very dense star. In Section IV, we present our conclusions. 



II. GINZBURG-LANDAU FREE ENERGY AND THE OPS 



In the superconducting phase the dominant pairing channel consists of two quarks of the same 
helicity. We will denote the corresponding order parameters by $l and $>r. &l,r are 3x3 matrices 
transforming by the 3 representation of SU(3)c and SU (3) One can argue that a semi-superfluid 
string has the same winding number for both and In our calculation we assume that 

&l = ®r = The Ginzburg-Landau(GL) free energy is a function of <E>. In the weak coupling and 
in the chiral limit it is given by [15| |9| 

T = \g%G% + ^F em F em + 2 KT tr(r3$) t (C$)+atr($ t $)+/3 1 (tr$ t $) 2 +/5 2 tr($ t $) 2 (1) 

If the static electromagnetic gauge field is A EM and the static color gauge fields are A a , a = 1, 8, 
the covariant derivatives and field strengths are 

l3$ = V$ - i g A a T a <S> - ieA EM T EM <f>, 
G% = d^-djAI+gf^A^ 
F EM = d % Af M -d Af M . (2) 

Now if the field is <& = $l(&r) then the corresponding flavor symmetry group is SU(3) F = 
SU(3) L (SU(3) R ). Under the element {V c ,V F ,e ioiB } € SU(3) C X SU(3) F x U(1) B , the diquark 
condensate transforms as 

$ -> V F <S>V$e laB (3) 

where Vf € (the 3 represenation of)SU(3)F, Vc € (the 3 represenation of)SU(3)c and e laB £ 
U(1)b- However from Eq.Q, we see that the group of elements {(zi, z^ 1 z^ 1 , z 2 ) : zi,z 2 — 
a cube root of unity} = Z3 x Z3 leaves $ invariant. So the symmetry group of the free energy 
r is 

^_ SU(3) c xSU(3) F xU(l) B 

G - z^Tz z ■ (4) 

In the symmetric or QGP phase, the diquark condensate vanishes and so is invariant under the 
group G. On the other hand $ is non-zero in the superconducting phase. As a result a smaller 
group {H G G) of transformations keeps $ invariant. The symmetry group H depends on the form 
of $ and thereby on the state of the superconducting phase. 

In the CFL phase the free energy T is minimized when $ is proportional to a constant unitary 
matrix U. So one can write for the minimum energy configuration, $0, 



$0 = r\U 



(5) 
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with rj is a positive real number. For the analysis of the symmetry breaking pattern one can take 
$o = r ?l without loss of generality, where 1 is a 3 x 3 identity matrix. The group ££7(3) x Z 3 : 
{(V, V z , z) : z = a cube root of unity, V G ££7(3)} keeps $o = tyl invariant. This set contains 
all the elements of Z 3 x Z 3 defined above. In order to find the stability group H C G of <J> we must 
quotient this set by Z3 x Z3. Hence 

H= sj M *z 3 

Z 3 xZ 3 y > 

The symmetry breaking pattern in the transition from normal to CFL phase (G — > H ) therefore is 

££7(3) c x ££7(3) F x U(1) B SU(3) x Z 3 



Z 3 xZ 3 Z 3 x Z 3 

Thus the order parameter space OPS for the <E> is given by, 



(7) 



OPS = [££7(3) c x ££7(3) F x £7(l) B ]/[££7(3) x Z 3 \ = £7(3) = SU ^* . ( 8 ) 

^3 

The OPS = £7(3) allows NNL's. In the language of homotopy groups, the NNL's are classified by 
the first homotopy group iri(OPS). In this case, 

iri(OPS) = Z. (9) 

We now explain its features. 

U(3) allows for non-abelian vortices. We can see this as follows. A nontrivial closed loop in 
U(3) can be described by a curve in SU(3) x U(l) beginning at identity which becomes closed on 
quotienting by Z 3 . Consider the curve from (1, 1) to (e* 2 ^ 3 , e^ 2 ^) in SU(3) x [7(1). In the SU(3) 
part, only the end points of this curve matter to determine the homotopy class of this curve. In the 
£7(1) part, it is the curve {e %v : < <p < — 2ir/3} from 1 to e~ i27r//3 in the anticlockwise direction. 
The S'£7(3) part of this curve is nontrivial. In £7(3) then, it is a nontrivial non-abelian closed loop /. 
It is the generator of Z in Eq.©. If [I] is the homotopy class of this loop, [I] 3 is associated with a 
closed loop in ££7(3) and £7(1). The closed loop in ££7(3) can be deformed to a point. Hence [I] 3 is 
a NNL in £7(1) and corresponds to the abelian superfluid vortices studied in ref.[8|- The elementary 
non-abelian vortices can be associated with [I] or Any non-abelian vortex is associated with 

[lW k or [l} 2 [l] 3k for k £ Z. 

Now we construct the loops in the OPS. We consider two loops, one in the homotopy class of [I] 
and the other in the homotopy class [l 1 ] = [l]~ 2 . The projections of these loops in the ££7(3) part 
of the OPS are same. The projections of loops I and I' in £7(1)b go from identity to e~ 27Tl ^ 3 in the 
anti-clockwise and clock-wise directions respectively. The latter explicitly is {e lip : < <p < 47r/3}. 
In U(1)b they are generated by the baryon charge, 

Qb = \l. (10) 
The loops I and I are parameterized in £7(1) b as 

e -iaQ B /2 e iuQ B /-q) 

where the parameter a varies from to 2tt. Both these loops start from identity and end at e 4 ™/ 3 
in£7(l) B . 
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Now let us discuss the projection of the loops I and I on the individual groups SU (3)c and 
SU(3) F . Note that the U{l) EM is a subgroup of SU(3) F . The generator T EM of U(1) BM in the 
(ds, su, it<i) basis is given by 




rpEM _ 

' 6 \ 1 J 

(12) 

which is linear combination of the generators S 3 and S 8 of SU(3) F . So a curve generated by T EM 
lies entirely in SU(3) F . To simplify the matter we work in a basis of generators such that T 8 (S 8 ) 
of SU(3) C (SU(3) F ) have the same matrix representation 




(13) 

as in ref.0. A path from (1,1) to e 27ri/3 (l,l) in SU(3) C X SL/(3) F can be generated by T EM 
or T 8 a linear combination of them. The projection of this curve on SU(3)c and on SU(3) F can 
be different. For example the projection of this curve can go from 1 to e 27 ™'' 3 in SU(3)c and be 
just a point in SU(3) F or vice- versa. Topologically there is no difference between these different 
possibilities. However dynamically they are different. When the projection in SU(3)c is from 1 to 
e 27 ™/ 3 and just a point in SU(3) F the resulting string configuration will be made of only color fields 
V3 - 

with g-^-<£> 8 = 27r. On the other hand if the projection in SU(3) F is from 1 to e 27Tl ^ 3 and is just 

a point in SU(3)c, the resulting string configuration will be made of only ordinary magnetic field 
with e$> EM = 2tt. But in CFL there is a mixing between A EM and A 8 into the following new gauge 
fields HU, 

A x = cos(A EM +sm(A 8 
Aq = - sin (A EM + cos C^ 8 (14) 

where Ax is massive and Aq is massless. C depends on the couplings as follows: 

C0S C = l/^ — 6 o i a ■ (15) 
y e 2 + 3g 2 /4 v ' 

Because of the mixing between the gauge fields a path from (1, 1) to e 27rl/3 (l, 1) in SU(3) C X SU(3) F 
has projection both in SU(3) F and SU (3)c for the minimum energy string configuration. As a result 
for minimum energy, the sum of the ordinary and color magnetic flux should satisfy 

(e$ EM )+g^ 8 = gA X =±2tt (16) 



where $ x is the magnetic flux of the massive gauge field A x and g x = y/ e 2 + 3g 2 /4. Note that even 
though the minimum energy configuration has a flux both of ordinary and color magnetic fields, a 
configuration with only ordinary magnetic flux can still be stable because of flux conservation. 

We mention here that previous studies 01 have considered the flux tube of the A x field. This 
flux tube is analogous to the electroweak string. Unlike our case these solutions correspond to a 
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closed loop in the SU(3) part of the OPS and are not topologically stable. Their counterparts in 
the electroweak theory have been extensively studied for stability. The results show that when 
the weak gauge coupling is larger than the abelian gauge coupling the solutions are unstable against 
expanding of the core of the string. For the same reason the flux tube [lfj considered previously in 
CFL will be unstable because the strong coupling constant is an order of magnitude larger than the 
electromagnetic coupling. 

We expect that in the gCFL phase also there will be semi-superfluid strings since the symmetry 
breaking pattern and the OPS are same as that of CFL case. In the mCFL phase we expect that the 
loops considered in the above discussion remain non-trivial. Since different components of $ have 
different mass in this case the core structure of the defect will be different from the case when quark 
masses are degenerate and are zero. The abelian loops [l] 3k as well as non-abelian loops [^][/] 3fc and 
[Z] 2 [i] still remain when quark masses are finite. Because of spontaneous breaking of U(1)b the 
OPS always contains a U(l). The effect of non-zero masses only changes part of the OPS which is 
generated by the nonabelian generators. However this change does not affect the above non-abelian 
NNS's as long as the nonabelian generator of this loops T a is not explicitly broken. So there should 
be non-abelian strings in the mCFL phase. The same arguments can be made about the dSC and 
uSC phase. For the 2SC and g2SC cases the situation is similar to the electroweak case. As the 
condensate has only one nonvanishing component, the loops generated by Qb and T 8 are same and 
hence will lie in SU(3)c,f- So one will not have any stable topological string solution in the 2SC 
and g2SC phases. 

III. FIELD EQUATIONS AND NUMERICAL SOLUTIONS FOR NA 
SEMI-SUPERFLUID STRINGS 

In this section we consider the I and I in the OPS and derive the action and the field equation for 
the semi-superfluid strings. As we mentioned above, because of the mixing between A EM and Ag, 
the NNL's will have projections in both SU(3)f and SU(3)c- In the following, we will denote the 
generator corresponding to the massive gauge field Ax by A. To simplify the notations we denote 
T EM by T in the following. The loops I and I' are parameterized by 



with a constant r/. 

The projection of {Mi (a) : a e [0, In]} in SU(3) goes from / to the center element e 27 ™/ 3 while 
the projection in U{1)b goes from I to e~ 27 ™/ 3 . The second loop {M 2 (a) : a 6 [0, 2tt}} is the same 
as Mi (a) except the path in U(1)b is covered with the reversed orientation. 

The minimum energy configuration strings are associated with Mi (a) or M2(a). To keep the 
energy contribution from the covariant derivative to its minimum, only the massive gauge field A is 
excited. Assuming that the string is along the z-axis and is cylindrically symmetrical around it, the 
$ configuration for the string (corresponding to Mi (a) with a = 8, 8 being the polar angle of the 
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position vector on the xy-plane) is given by 

/ vf(r)e' me \ [ 4>(r,6) \ 
$(r,0) = ^i(r) = ViW 

V ^,(r) / V Mr) J 

(17) 

To simplify our calculations we assume that ipx%{r) = r/. The string configuration with the lowest 
energy will have some r dependent profile for ipi.2 (r). Near the core of the string ?/>i.2 will have 
values slightly larger than r\ due to coupling with the <fr field and the gauge fields . Far from the 
core of the string where r — > 00, ipi^ir) must be equal to i] but without any nontrivial winding like 
the <f> field. 

The finiteness of the potential energy part of the free energy (Eq.Q) requires that f(r — ► 00) = 1. 
The kinetic energy part of the free energy can be minimized by an appropriate choice of the gauge 
fields A. Since the string is cylindrically symmetrical, the phase varies only in the 9 direction in the 
x, y plane, so only the Ag will be non-zero. The total static free energy corresponding to the loop 
{Mi (a), a £ [0, 2tt]} with cf> = r}f(r)e~ ine is given by 



r(4>, A) = 2n T \(d + l ^A)^ + ^E± A ^ + - d y A,f 

+a(|0| 2 + 2t7 2 ) + /3(2</> 4 + Aifcj? + 6r? 4 ). 



(18) 



The above free energy is minimized if, at a large distance r from the string, Ag is a pure gauge field. 
This implies the following form for the gauge field Ag: 



Ag = ^-,j(r^ 00) = 1. 



(19) 



Now the covariant derivative in the 9 direction is given by 



D e ^ = d e $ - ig x A e T<S> = 77 







/ .n + 2ig^ e Q 
r 3r 









3r 




3r 



(20) 



The gradient energy density from the variation of the fields along the 9 direction is given by 

" [(3n - 2g x tf + 2(g x £) 



\Dg<£\ 2 - 77 <> J " • " J 



9r 2 

which is minimized by g x t; = 1 for n — 1. So the gauge field Ag takes the form 



(21) 



Ag 



9xT 



(22) 



It is important to note that the gauge field reduces the gradient energy by — for the string 

o 

configuration corresponding to loop M\ (9) . 
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The Euler-Lagrange equations for <I> and the gauge field A respectively are 



o A,. + — - — (0 a M - ) „ — !</>! ^ „ ^ M = 0.(24) 



The static equations satisfied by / and 7 are 



/"(r) + r/'W- 



f(r) ( 4 7 (r) 



/2/ 2 (r) + 1 



-53 — /( r ) -—vf ( r ) = 0, 

8 KT3 ^V 2 (r) 



7 (r) - -7 (r) - 8 KTff 2 7? 2 ^ - v 'g ' ^ ) 7(0 + — — 1 = 



(25) 
(26) 



Now let us consider the string solution for the loop given by M 2 {9). The ansatz for the $ field in 
this case is 



ip(r) 
$(r,0) = I r]f(r)e 



inO 






V f(r)e l 



(27) 



As in the previous case of the string solution corresponding to the loop Mi, we assume that 
ip(r) = 77 to simplify our calculations. The free energy with <fi = rjf(r)e m is given by 



r($, A) = 4k t \0- ^fAW + + J _ d v An) 2 

+a(2\cj)\ 2 + if) + /3(60 4 + 4r/ 2 2 + 2T? 4 ). 



(28) 



Again since the phase of the condensate depends only on 9, only Ag is non-zero. To find out the 
asymptotic form of Ag, we impose the condition that it asymptotically becomes a pure gauge field 
and also minimizes the gradient energy. The covariant derivative of $ is 



D e <f> = dg§ - ig x A B T§ = rj 



( 2»gs£ 
3r 








, n _ igJi_ in0 







{j n ig^ g 
r 3r ' I 



The gradient energy density corresponding to the phase variation of $ is now given by 



(29) 



\D B <S>\ 2 = — [2(Zn- g x Z) 2 + A{ 9x 2 ] 



(30) 



which is again minimized by g x £ = 1. The reduction in the gradient energy due to the gauge fields 

2 

is now only by a factor of - . The Euler-Lagrange equations for the fields 4> and A are 
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d »mm ~ w = 4kt(5m " l ^ A ^ ~ l -f A ^ - 2(s + 2 ^ - 1 W|2< ^ = °« (31) 

d "WA^-SA^ = 3— (^^-W)--g-^ ^-—9-^ = (32) 

which simplifies to the following: 



r r"* \ 3 / 2kt kt 
7 W - -7 M - 8K T g K ?7 i l(r) + = 0. (34) 



(33) 



We solve equations J5^J-J5BJ| and l|3"3"|) - (|3"P numerically to find out the string profile for the two 
loops I and I' respectively. We choose the values of parameters which will correspond to a type II 
color superconductor. For our calculations we took g x = 2.0, kt — 0.42, — 1.26 and rj = lOOMeV. 
The parameter a is obtained from the relation a — SPrj 2 = 1.008 x 10 5 MeV 2 . We also use the 
following relations 0: 

a (3ri 2 2 
m 

2K T glrf = m\. 



(35) 



For these parameters we have the Higgs mass m$ ~ 245MeV and the Meissner mass (inverse of the 
penetration depth) vtla ~ 183MeV. 

The results of our numerical solutions are shown in Fig.l and Fig. 2. Though the figures look 
similar, the energy of the configurations are very different. One can see that both the <fi and A 
profiles vary more slowly and reach their asymptotic values at larger r for the string corresponding 
to the NNL I'. The NNL I' in the OPS travels a longer path in U(1)b which costs larger gradient 
energy for the semi-superfluid solutions as U(l) b is a global symmetry group. 

In the profiles of the string, we see that far from the core of the semi-superfluid string, the field 

is color-flavor locked, i.e. the field is a constant times a SU (3) matrix. However at the core of the 

string the condensate is not locked. 

The NA string corresponding to the loop Mi (a) 3 and a winding number one superfluid string are 

3 

topologically equivalent. For both these string configurations the energy density is ~ at large 

r. It is not clear if one of these configuration decays into the other. On the other hand, this NA 
string is also topologically equivalent to three elementary NA strings. However we do not believe 
that three clearly separated NA elementary strings will evolve into the above mentioned NA string 
or the superfluid string. We expect that elementary superfluid strings of the same winding number 
repel each other at large separations. So the total flux of a network of clearly separated strings 
should be additive. We address this issue in a future work. 



IV. FORMATION AND IMPLICATION OF STRINGS IN CFL PHASE 



Our non-abelian string solutions are the only stable topological string solutions with quantized 
flux of ordinary as well as color magnetic field. All other solutions with flux tube are topologically 
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FIG. 1: The /(r) and 7(7") profiles for the string corresponding to the non-trivial loop I. 



unstable for realistic values of the strong and electromagnetic couplings. It is interesting to note that 
when one takes an electron in a closed path around the non-abelian strings, the Aharonov-Bohm 
phase is different from 2tt. This will lead to strong scattering of electrons from these strings. In the 
following we discuss a possible scenario for the formation of string defects inside the core of very 
dense stars. 

Semi-superfluid or superfluid defects in the CFL phase can form during the phase transition from 
QGP to CFL phase. Topological strings can also be induced in the CFL phase from the outer 
surrounding confining medium when the star starts to spin up. This is in analogy with creation 
of vortices in rotating superfluid. In the following, we discuss the formation of strings during the 
QGP-CFL transition inside very dense stars. 

A phase transition from the QGP to CFL phase may be expected inside the core of dense stars 
while the star is cooling. Topological strings as well as non-topological strings will be formed during 
this transition. We assume that this transition is of first order. In this case, the transition takes 
place via nucleation of CFL bubbles in the QGP background. When the temperature inside the star 
cools below the critical temperature T c , bubbles of CFL phase nucleate in the QGP background. 
Inside the bubbles, the magnitude of the condensate is uniform and the massive gauge field A is 
zero. The phase a of the condensate is also uniform inside the bubble and varies randomly from one 
bubble to the other. When the bubbles meet, a interpolates between the values of the phase in the 
two colliding bubbles. There are two possible ways a can interpolate in the presence of gauge fields 
[20| . Numerical ex peri ments have shown that the interpolation of a is such that the total variation 
of a is the lowest [2l[. When a starts to interpolate, the massive gauge field A also gets excited 
to minimize the gradient energy. When three or more bubbles collide, it sometimes leads to the 
variation of a by an integer multiple of 2n around the loop at the intersection point. When this 
happens, a semi-superfluid string is formed. This is the conventional mechanism of defect formation 
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FIG. 2: f(r) and j(r) profile for string corresponding to the non-trivial loop l' 

known as the Kibble mechanism |Tflj . However there are other mechanisms which also contribute to 
the defect density [22|. 

One may expect that the presence of strong external magnetic fields will affect the formation of 
semi-superfluid strings. However since the unbroken gauge field in the CFL phase consists ~ 99% 
of the electromagnetic field, the major part of the external magnetic field will propagate unscreened 
[23| . Only a very small fraction of the external magnetic field, basically related to the massive gauge 
field, will be repelled by the CFL phase. So in a sense the formation of semi-superfluid strings is 
more like the formation of flux tubes in ordinary superconductors under a small external field. Note 
that the massive gauge field is made up mostly of the color gauge field, so our flux tube will consist 
mainly of color magnetic flux. We expect that there is no long range color magnetic field in the 
medium before the transition takes place. So the formation of semi-superfluid strings during the 
transition is spontaneously induced rather than external field-induced. 

Now we discuss a possible scenario for the network of strings inside the CFL core of the dense star. 
Usually the temperature at the center of the star is higher and gradually decreases as one moves 
radially outwards. So when the star cools, the transition from the quark-gluon plasma (QGP) to 
the CFL phase will take place first in a thin spherical shell where the temperature drops below 
the critical temperature. In this thin spherical region bubbles of CFL nucleate and grow. As they 
grow the bubbles coalesce with other bubbles. However since the temperature is higher towards the 
center of the star the bubbles will grow mostly in the spherical region forming a thin spherical shell 
of CFL. The picture of the phase distribution is that of a spherical shell of the CFL phase covering 
the QGP core with temperature above T c . The QGP and CFL phases are separated by the QGP- 
CFL boundary. The CFL shell and the outer confined crust are separated by the CFL-confining 
boundary. At this point, the cooling of the star will be different from that of cooling due to neutrino 
emission because there will be generation of latent heat when QGP converts into CFL. 

Further dynamics of the transition can be cither by bubble nucleation or motion of the interior 
wall of the spherical CFL shell covering QGP. However these two pictures are not much different 
because even if there is bubble nucleation, the bubbles will nucleate close to the boundary wall as 
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the temperature farther inside is either close to T c or higher. So the basic picture of transition is that 
a CFL shell is formed due to bubble nucleation and then the phase transition takes place through 
the motion of the interior wall of the CFL shell towards the center of the star. 

The strings which can survive such a transition are those which are oriented along the radial 
direction. One end of this string will end on the inner QGP-CFL boundary and the other in 
the CFL-confining boundary. The strings can end on the CFL-confining boundary because of the 
availability of color monopole and anti-monopoles pairs in the confined phase . The strings with 
both ends connected to the outer confining crust will decay by shrinking to the crust. The surviving 
radially oriented strings will possibly increase in length along the radial direction. However the 
density of string ends on the inner CFL-QGP boundary will increase due to the shrinking QGP 
core. Eventually the density will be so high that ends of different strings will come in contact with 
each other. One can argue that the total number of strings ending on the QGP-CFL boundary 
are even in number with equal number of strings and anti-strings. Two such strings with opposite 
winding will join together forming huge ^-shaped strings with ends connected to the outer confining 
crust. These ^-shaped strings are unstable to moving towards the outer confining crust. This 
movement happens almost simultaneously to all string-antistrings pairs ending on the QGP-CFL 
boundary. Such an evolution of a network of strings may affect the properties of a star like its 
angular momentum. 

V. CONCLUSIONS 

We have studied non-abelian semi-superfluid strings in high density QCD. Inclusion of gauge fields 
reduces the energy of these strings compared with the U(1)b superfluid string. Even with the gauge 
fields the energy per unit length of the semi-superfluid string is logarithmically divergent with the 
system size. Still such strings are relevant for finite system like stars. This is unlike the flux tubes in 
ordinary superconductors. The semi-superfluid strings are partly like superfluid strings and partly 
like flux tubes. These are the only topological strings possible in high density QCD which have flux 
of ordinary magnetic and color magnetic fields. Parallel transport of an electron in a closed loop 
around these strings picks up an Aharonov-Bohm phase different from 2ir leading to their strong 
scattering from the semi-superfluid string. We propose a scenario for a string network at the core 
of the star, the evolution of which can affect the dynamics of the star. However a detailed study 
using realistic phase transition dynamics for the QGP-CFL transition and formation of strings is 
necessary to make any definite prediction for the string network or their possible effects. We plan 
to do such calculations in the future. 
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